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Abstract. We review the determination of the NLO Brodsky-Lepage-Mackenzie (BLM) 
renormalization scale for the pion transition form factor. We argue that the prediction 
for the pion transition form factor is independent of the factorization scale at every 
order in the strong couphng constant. 

1 Introduction 

The pion transition form factor, the simplest exclusive quantity, offers an excel- 
lent testing ground for QCD. For large virtualities of the photons (or at least 
for one of them) perturbative QCD (PQCD) is applicable A specific feature 
of this process is that the leading-order (LO) prediction is zeroth order in the 
QCD coupling constant, and one expects that PQCD for this process may work 
at accessible values of spacelike photon virtualities. 

The pion transition form factor F^.^(Cp) is defined in terms of the amplitude 

r^^- ^le^ F^^{Q^)e^'''''^ P^qp, (1) 

and for large-momentum transfer — ~q^, it can be represented j|,|l| as a 
convolution 

F^^{Q'')^'P*[x,^il) (g, TH{x,Q^f4), (2) 

where (g) stands for the usual convolution symbol {A{z)i^B{z) — dzA{z)B{z) ). 
In (|^), the function Th{x, Q^, fJ-^p) is the hard-scattering amplitude for producing 
a collinear qq pair from the initial photon pair; 'P*{x, fi'^p) is the pion distribution 
amplitude (DA) representing the probability amplitude for finding the valence qq 
Fock state in the final pion with the constituents carrying fractions x and (1 — x) 
of the meson's total momentum P; fi^p is the factorization (or separation) scale 
at which soft and hard physics factorize. In this standard hard-scattering ap- 
proach, pion is regarded as consisting only of valence Fock states, transverse 
quark momenta are neglected as well as quark masses. 
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The hard-scattering amphtude Tu is obtained by evahiating the 7*7 — > qq 
amplitude, and has a well-defined expansion in Q!s(/i|j), with being the renor- 
malization (or coupling constant) scale of the hard-scattering amplitude. Thus, 
one can write 

+ ^^Tl^\x,Q',^^%^^,l) + ■■■. (3) 

The process-independent function <P(x,^^) is intrinsically nonperturbative, 
but it satisfies an evolution equation of the form 

d 

where V{x,u,as{lJ.%)) is the perturbatively calculable evolution kernel. If the 
distribution amplitude <?(x,/Kq) is determined at an initial momentum scale 
(using some nonperturbative methods), then the differential-integral evolution 
equation (^) can be integrated using the moment method to give <P{x,fi^). 
The perturbative expansion of the pion transition form factor takes the form 

F,.{Q') = F\l\Q') + fW(Q^) + F(lHQ^^.l) + • • • . (5) 

The choice of the expansion parameter represents the major ambiguity in 
the interpretation of the perturbative QCD predictions. We see that the cou- 
pling constant as{n%), as well as, the coefficients Fj^ [i > 1) from (||), depend 
on the definition of the renormalization scale and scheme. The truncation of 
the perturbative expansion at any finite order causes the residual dependence 
of the prediction on the choice of the renormalization scale and scheme, and 
introduces the theoretical uncertainty. If one can optimize the choices of the 
scale and scheme according to some sensible criteria, the size of the higher-order 
correction as well as the size of the expansion parameter, i.e. the QCD running 
coupling constant, can then serve as sensible indicators of the convergence of the 
perturbative expansion. 

The simplest and widely used choice (the justification for the use of which 
is mainly pragmatic), is to take the scale equal to characteristic momentum 
transfer of the process, i.e. in our case = Q^. But since each external mo- 
mentum entering an exclusive reaction is partitioned among many propagators of 
the underlying hard-scattering amplitude, the physical scales that control these 
processes are inevitably much softer than the overall momentum transfer. 

Several scale setting procedure were proposed in the literature H:^,^ . In the 
Brodsky-Lepage-Mackenzie (BLM) procedure js) , all vacuum-polarization effects 
from the QCD /3- function are resummed into the running coupling constant. Ac- 
cording to BLM procedure, the renormalization scale best suited to a particular 
process in a given order can be, in practice, determined by computing vacuum- 
polarization insertions in the diagrams of that order, and by setting the scale 
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Fig. 1. Pictorial representation of the pion transition form factor calculational ingre- 
dients. 

demanding that n/-proportional terms should vanish. The optimization of the 
renormahzation scale and scheme for exclusive processes by employing the BLM 
scale fixing was elaborated in Q . The renormahzation scales in the BLM method 
are physical in the sense that they reflect the mean virtuality of the gluon prop- 
agators and the important advantage of this method is "pre-summing" the large 
(Pocts)^ terms, i.e., the infrared renormalons associated with coupling constant 
renormahzation (|^ and references therein). 

In our recent work we have determined the BLM scale for the pion tran- 
sition form factor, i.e., for the 7*7 — > tt process. The LO prediction for the 
pion transition form factor is zeroth order in the QCD coupling constant, the 
NLO corrections Q represent leading QCD corrections and the vacuum polariza- 
tion contributions appearing at the next-to-next-to-leading order (NNLO) were 
needed to fix the BLM scale from from the requirement 

J^^''"'^(Q',Mlf, = MlLM)=0, (6) 

where -F7^'"^''(Q^, represents the ri/-proportional NNLO term from (||). 

In this work we outline important points of this calculation and present the 
results that follow from the consistent calculation up to n/-proportional NNLO 
contributions to both the hard-scattering and the distribution amplitude. 

2 Analytical calculation 

We first outline the calculational procedure and its ingredients which are illus- 
trated in Fig. ^ 

The 7* + 7 — 99 amplitude denoted by T contains coUinear singularities, 
and it factorizes as 

T{u,Q^)=TH{x,Q\fil) ® ZTMx,u;fil). (7) 
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Here, fij^ denotes a factorization scale at which the separation of coUinear singu- 
larities takes place, and all coUinear singularities are factorized in Zt^coI, since 
Th is, by definition, a finite quantity. 

On the other hand, a process-independent distribution amplitude for a pion 
in a frame where P+ = P'^ + = 1, p- = P° - = 0, and Pl = is defined 
as 



(8) 

where i? = exp is a path-ordered factor making ^ gauge 

invariant. The unrenormalized pion distribution amplitude <P{u) given in (^) and 
the distribution amplitude renormalized at the scale are related by 

a multiplicative renormalizability equation 

^{u)=Z^^ren{u,V-H%)(^^{v,nl,). (9) 

By convoluting the "unrenormalized" (in the sense of coUinear singularities) 
hard-scattering amplitude r(w, Q^) with the unrenormalized pion distribution 
amplitude 'P{u), given by (0) and (P), respectively, one obtains 

F^^(Q2) ^ <pt(y) ^ T{u, Q2) . (10) 
The divergences of T{u, Q^) and ^(u) cancel 

Zt,coi{x, U] flp) (g) Z^^reniu, V\ ^ip) = (5(x - u) , (11) 

and the usual expression (|^) emerges. It is worth pointing out that the scale fj,"^ 
representing the boundary between the low- and high-energy parts in (|^) is, at 
same time, the separation scale for coUinear singularities in T{u, Q^), on the one 
hand, and the renormalization scale for UV singularities appearing in <P{u), on 
the other hand. 

We note also that the pion distribution amplitude as given in (|8|), with |7r) 
being the physical pion state, of course, cannot be determined using perturbation 
theory. We can write ^(u) as 

<P{u) ^ ^{u,t) ® (qq;t\TT) , (12) 

where (f>{u,t) is obtained from (||) by replacing the meson state \n) by a \qq]t) 
state composed of a free quark and antiquark. The amplitude (f) can be treated 
perturbatively, making it possible to investigate the high-energy tail of the pion 
DA, to obtain Z^^^en and to determine the DA evolution. 

We proceed to calculation. This is the first calculation of the hard-scattering 
amplitude T(u, Q^) of an exclusive process with the NNLO terms taken into 
account. The subtraction (separation) of coUinear divergences at the NNLO is 
significantly more demanding than that at the NLO. Owing to the fact that the 
process under consideration contains one pseudoscalar meson, the calculation is 
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further complicated by the 75 ambiguity related to the use of the dimensional 
regularization method to treat UV and coUinear divergences. The consistent 
calculation of T and 4> enable us to resolve these problems and, hence, we have 
calculated the LO, NLO, and n /-proportional NNLO contributions to the pertur- 
bative expansions of both the hard-scattering amplitude and the perturbatively 
calculable part of the distribution amplitude. 



3 Discussing the factorization scale independence of the 
finite order result 

The dependence of pion distribution amplitude ^{Xji^jp) on fij^ is specified by 
the evolution equation (Q). This dependence is completely contained in the evo- 
lutional part (pv 

${v, ^1%) = (t)v{v, s; 11%, III) ${s, III) , (13) 
which satisfies the evolutional equation 



/J-F-^^t^viv^S^lip^Ill) = V{v,s\lJ.%) (g) 0y(s',S,/i|,^o) , (14) 
Ullp 

while <?(s,/i§) represents the nonperturbative input determined at the scale /Xq. 

By differentiating (^ with respect to and by taking into account (Q), one 
finds that the hard-scattering amplitude satisfies the evolution equation 

tJ-l^THix,Q'\iil)^~THiy.Q^li%) ® V{y,x;iJ.%), (15) 

which is similar to The dependence of Th{x, Q^,ii%) can be, analogous 
to (|l3|), factorized in the function (pv (y , x , Q'^ , n\) as 

TH{x,Q\iJil)^TH{y,Q\lil^Q'')®4>v{y.x,Q\t^l)- (16) 

Using ( p^ one can show by partial integration that (|l^) indeed represents the 
solution of the evolution equation (Esh. 

By substituting ( ^ and (|l6| ) iii(||), we obtain 

F^,(Q2) = Th{v, g^ Q') ® 0y(y, s, Q\ 111) ® <Z'*(s, i^l) , (17) 

where 

<Pv{y, X, Q^, lip) «i (f>v{x, s, lip, III) = (j)v{y, s, 111) , (18) 
has been taken into account, ft is important to realize that the expression ( p^ 
is valid at every order of a PQCD calculation, and this can be easily shown (see 
0). Hence, the factorization scale disappears from the final prediction at 
every order in 0^5 and therefore does not introduce any theoretical uncertainty. 
The crucial point is that both the resummation of {as \n{ii'^ / nl))" terms in <P as 
well as the resummation of {as ln((5^//x|,))" terms in Th, have to be performed 
using ( [1^ ) and ( |l^ ) along with the results from (14). We note here that by 



adopting the common choice 11^ = Q^, we avoid the need for the resummation 
of the {as ln((5^//i|,))" terms in the hard-scattering part, making the calculation 
simpler. 
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Fig. 2. The LO and NLO predictions for the pion transition form factor obtained using 
the MS scheme (and the usual one-loop formula for as with Ajrj = 0.2 GeV^). 



4 Numerical predictions 

We refer to for the complete analytical expressions for the pion transition 
form factor calculated up to n/ proportional NNLO terms. 

The prediction for the pion transition form factor and the BLM scale 
depend on the form of the distribution amplitude. There is increasing theoret- 
ical evidence coming from different calculations [ p^ that the low energy pion 
distribution amplitude does not differ much from its asymptotic form. 

The expression for the pion transition form factor Q"^ F^t^{Q'^) corresponding 
to the asymptotic distribution reads 



47r (47r)2 

2 \ / „2 
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-43.47- 20 In 



(19) 



where = y/2/(!> is a flavour factor, while /^r = 0.131 GeV. The n/-proportional 
NNLO contribution determines the value of the BLM scale 

^^'n-{ALMr^^■ (20) 

One notes that this scale is considerably softer than the total momentum transfer 
Q^, which is consistent with partitioning of among the pion constituents. 

The NLO predictions obtained in the MS scheme are displayed in Fig. ^. The 
predictions based on the asymptotic DA are, in contrast to the ones obtained 
using the CZ DA in good agreement with the experimental data |p^ . 

Nevertheless, the rather low BLM scale given in (^0|), and consequently the 
large ^^(/^bla/)' questions the applicability of the perturbative prediction at 
experimentally accessible momentum transfers. The NLO predictions obtained 



BLM scale for the pion transition form factor 



7 



0.5 



0.4 



0.3 

O 



■ CLEO(1998) 



S 0.2 



0.1 



0.0 



8 12 

Q' [GeV'l 



16 



20 



Fig. 3. The LO and NLO predictions for the pion transition form factor in the av 
scheme {Av = 0.16 GeV^). 



assuming the asymptotic DA and the BLM scale ( po|) satisfy the requirement 

"s(/^fl) < 0.5 for Q2 > 6 GeV^. This rehability hmit is indicated on Fig. | The 

transition to the more physical ay scheme, may offer a way out of this problem. 
In the exclusive hadronic amplitudes were analysed in the ay scheme, in 

which the effective coupling ay(/i^) is defined from the heavy-quark potential 
The av scheme is a natural, physically based scheme, which by definition 

automatically incorporates vacuum polarization effects into the coupling. The /iy 

scale reflects the mean virtuality of the exchanged gluons. 

If use is made of the scale-fixed relation between the couplings otj^^ and cty 
then, to the order we are calculating, the NLO prediction in the av scheme 

is obtained by taking = /iy = e^/"^ I^%laI' i-*^- f'^'' ^^'^ asymptotic DA 

{^^'vr^e'^H^^lLMr^^■ (21) 

The NLO prediction for Q^FyTy{Q^) obtained in av scheme is depicted in Fig. |[ 
As can be seen, it is in good agreement with experimental data. We note that, 
since av is an effective running coupling defined from the physical observable, 
it must be finite at low momenta, and the appropriate parameterization of the 
low-energy region should in principle be included (see [ p^ for various proposals). 

5 Conclusions 

In this paper we have reviewed the determination of the NLO BLM scale for the 
pion transition form factor. A consistent calculation of both the hard-scattering 
and the perturbatively calculable part of the pion distribution amplitude has 
been performed up to n/-proportional NNLO terms. 

It has been demonstrated that the prediction for the pion transition form 
factor is independent of the factorization scale at every order in the strong 
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coupling constant as- Provided both the hard-scattering and the distribution 
amphtude are treated consistently regarding their dependence, the factor- 
ization scale disappears from the final prediction at every order in as without 
introducing any theoretical uncertainty. One can use = to simplify the 
calculation, but any other choice would lead to the same result. 

The renormalization scale nj^ fixed according to the BLM scale setting pre- 
scription within the MS scheme and corresponding to the asymptotic pion distri- 
bution amplitude, turns out to be Mbla/ ~ Thus, in the region of < 8 
GeV^, in which the experimental data exist, /i^LM < ^ GeV^. Consequently, 
the prediction obtained with = A^^LAf cannot, in this region, be considered 
reliable. 

In addition to the results calculated in the MS renormalization scheme, the 
numerical prediction assuming the same distribution but in the ay scheme, with 
the renormalization scale /i^ = jiy = e^^^fJ-%LM ~ has also been obtained. 

It is displayed in Fig. js] and, as seen, is in good agreement with experimental 
data. Due to the fact that the scale fiy reflects the mean gluon momentum in 
the NLO diagrams, it is to be expected that the higher-order QCD corrections 
are minimized, so that the leading order QCD term gives a good approximation 
to the complete sum. 
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